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Abstract 



We present results of a set of experiments to investigate the effect of dis- 
sipative external electromagnetic environment on tunneling in linear arrays 



o 

O ' of junctions in the weak tunneling regime. The influence of this resistance 



pq ' decreases as the number of junctions in the chain increases and ultimately 

>' 
OO ' becomes negligible. Further, there is a value of external impedance, typically 

cn ; 

^i^ I ~ 0.5 ki7, at which the half-width of the zero-voltage dip in the conductance 

Q\ ' curve shows a maximum. Some new analytical formulae, based on the phase- 

^' 

correlation theory, along with numerical results will be presented. 
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In recent years large attention has been paid to the role of electromagnetic environment 
on charging effects in small tunnel junctions, both theoretically and experimentally [|1|- [0. 
Yet, arrays of such tunnel junctions with well-defined external impedances have not been 
extensively discussed. This is partly because the theoretical formulation of such arrays 
is more elaborate: there are, e.g., cotunneling effects, inhomogeneities, and background 
charges, which along with the effects of environment make such an investigation rather 
difficult in general terms. This lack of theoretical predictions, in turn, has decelerated 
experimental search for observation of new features in arrays. In this letter we will attempt 
to fill part of this gap by demonstrating a set of experimental observations and a comparison 
of them to the results obtained from the already existing phase-correlation (PC) theory for 
single tunnel junctions, which we have now extended to analyze junction arrays numerically. 
This analysis is important in setting limits to the systematic error of the reading of the 
Coulomb blockade primary thermometer |^. 

According to the PC theory, the tunneling rate through the hh junction of a completely 
symmetric array with Ck = C, RT,k = Rt and Co,fc = |Q (see Fig. 1), in the weak- 
tunneling regime, Rt^u ^ Rk = h/e^, can be written as a convolution integral of the form 
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Vt{{n}) ^ V{5Ft {n}) = -^ / dE Pu{-5F^ - E). (1) 
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Here SF^~^ is the change in free energy of the array when an electron tunnels to right 
(left), {n} = {ni,n2, . . . ,njv_i} designates the charge configuration on the islands, and 
Pk{E) = {2txTi)~ j^^ dte'-'^^^'>^^n* is the probability density for the electron to exchange 
energy E with the environment. The correlation function Jk{t)i which accounts for the 
environment of the Wn junction, is given by 






where Z^ioj) is the total impedance of the circuit as seen by this junction, and (3 = {ksT) ^. 
By applying Fourier transform techniques to Eq. (1) one obtains 



r± tl 



This equation is central in the following numerical calculations. With F^'s and the algorithm 
in []TU| with i?s = Y^k=i^T,k-, "we find the current / through the array in equilibrium as 



X^i visited 
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In the expression above, starting initially from an arbitrary configuration, the states 
visited, {n}j,jsjjed, over which the outer sums run through, are obtained by dividing the 
interval [0, 1] into segments proportional to T^'s in each current state, and drawing a random 
number r in the interval. For sequential tunneling, the segment to which r corresponds to 
will specify the junction through which the tunneling event happens and the tunneling 
direction. This way the distribution of {n}msited will be statistically collected, and it will 
allow one to calculate the sums (now weighted by the distribution) according to Eq. (|[) 
similarly to what has been done in [|T0|. 

In the case of a symmetric two-junction array we will also use a simpler algorithm, 
described in [0], to obtain the probability of finding n excess electrons on the island, a{{n}), 
and finally the equilibrium current / through the array 

oo 

/ = 4 = e E ai{n})[rti{n}) - T, ({n})]. (5) 
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Assuming a completely symmetric array and a purely resistive environment Re, the com- 
mon real part of the total impedance, Zt{u!) = Z^{uj), in Eq. (2) reduces to the simple 
form Re[Zi(a;)] = Re/[{to/uJcY + N'^], with tUc = l/ReC. It is already suggested by this 
equation that the effect of the environment decreases with increasing A^ and becomes van- 
ishingly small for long arrays. Furthermore, using the partial sum expansion of coth(x) (or 
by applying Cauchy's integral theorem), one can evaluate J(t) = Jk(t) in Eq. (2) as 

M^ Rk 2 pn :^-^ I-Knlujn^ - (Awcj^J 



where uJn = 27m/ (3h are Matsubara frequencies. The above equation is a straightforward 
extension of the result obtained for a single tunnel junction in resistive environment P,p!T|. 
In our numerical calculations we have used the generally valid formula in Eq. (2), but 
equivalence of results is checked both by direct comparison of the numerical values derived 
from Eqs. (2) and (6), and by comparing the final conductance curves; the results from the 
two methods are indistinguishable. 

Figure 1 shows a schematic view of an array with its bias circuitry. To check for the 
consistency of results for arrays with different number of junctions, each sample included 
one pair of arrays with about 3 jum space in between. Each pair had a different number 
of junctions, typically A^ = 1, 2, A^ = 2, 8, and A^ = 2, 20 (only samples with A^ = 2, 8 
and A^ = 2, 20 are shown in Table I). The array consisted of Al/AlOx/Al tunnel junctions 
(0.01 — 0.05 fim^) with four chromium resistors, Z^j^uj) = Rj, at a distance of 2 /im, two 
at each end of the array. Cr resistances of 1 — 20 kQ can be easily obtained by adjusting 
the width (~ 100 nm), length (~ 2 fim), and thickness (~ 3 — 8 nm) of chromium films. 
The equivalent environment resistance {additional to the natural free space like impedance 
of ~ 100 Q) will, then, be Re = RiR2/{Ri + R2) + RsR^/iRs + -^4)- Samples were made by 
e-beam lithography and three-angle evaporation techniques. All measurements were carried 
out at 4.2 K. More details of measurement techniques are presented in 0. 

Figure 2 shows measured data for a typical sample with N = 2, Rt = 44.9 kQ, R^ = 1.12 
kQ and C = 2.25 fF (sample 9A in Table I (a)). C was determined by searching for the best 
fit of the measured depth of the zero-bias anomaly to the value given by Eqs. (1 — 3); those 
equations are in agreement with the zero-bias results of [Q. Each point in the simulated IV- 
curve (not shown in the figure) was obtained as a result of 1000 (pseudo) random tunneling 
events. Conductance curve (CC) was obtained by numerical derivation of the Jl^-curve. It 
is worth mentioning that even with much smaller number of draws, e.g. 10, the result agrees 
within better than 1.5% accuracy (with respect to the height and width of the CC dip) with 
those obtained from "long simulations". The time needed for each step, with a reasonable 
partition of voltage interval and using a regular desktop computer, was about 1 minute. 



The result with N = 2 obtained from Eq. (5), based on the algorithm presented in JTO|, is 
identical with that of a more comprehensive method described above. 

In Fig. 3 we have drawn the normalized half-width of the zero-bias minimum of the 
conductance curve, Vl/2, as a function of the environment resistance for different pairs of 
samples. By normalization we mean scaling the half-widths by those derived for arrays 
without an external impedance, i.e., scaling by Vi/2,0 = 5.4:39NkBT/e [§]. Usually, we have 
done measurements for each sample in six different combinations of current and voltage 
probes (four-probe measurement) and there are variations between the different combina- 
tions shown by the error bars. The unequal height of the error bars for different samples 
is due to this. For N = 2, V1/2 shows a sharp maximum at around Re — 1 kQ. For longer 
arrays, A^ = 8 and A^ = 20, the dependence is much weaker, and V1/2 stays close to V^i/2,0) 
which is directly demonstrating the advantage of using long arrays in Coulomb blockade 
thermometry. In each case (unambiquously only for N = 2, though), the normalized V1/2 
approaches unity with Re ^ and for R^ —>■ 00. This is in agreement with predictions of 
the PC theory. We will discuss this in further detail below. 

In the same figure the results obtained by numerical simulation, and for A^ = 2 by the 
direct calculation described above, are depicted. In spite of the overall agreement between 
experiment and theory, there is a noticeable discrepancy between the measured and the 
predicted value of Vi/2,max, the widest half- width, of the two-junction array. In this case 
(A^ = 2) we could calculate the half-width with the two methods described above, without a 
noticeable difference between them. For A^ = 20 the absolute difference between experiment 
and theory is smaller; the predicted peak itself is very small, less than 1% (see inset of Fig. 
3). Comparison of experimental data to those derived from the theory for N = 2 shows 
that the former produce even 15% wider conductance curves (Fig. 3). Cotunneling and 
other higher order tunneling effects may play some role in our samples, because the junction 
resistances are not large as compared to Rk- While in the simpler case of a two-junction 
array without electromagnetic environment higher order tunneling tends to broaden the 



half-width of the CC dip [0 , making such comparison between our data and the theoretical 



study of cotunneling effects in an array with dissipative environment [0 is beyond the scope 



of this work and is not done here. There is also a difference in the shape of the "shoulders" 
of the measured conductance curves for A^ = 2 in particular (see Fig. 2) as compared to 
the theory. Such a distortion of shape can be caused by the nonuniform size distribution of 



junctions in the array, which we have studied in the case of no external impedance P,P 
Experimentally the intercomparison of the depths in different samples is difficult because 
the size of the junctions varies from sample to sample, and this gives the main contribution 
to the depth variation. 

Next, let us consider the conductance curve in more detail. Using the time-domain 
formulation of single electron tunneling presented in PJT5[] we repeated calculation of the 



high-temperature conductance G{V) of the symmetric A^-junction array without stray ca- 
pacitance in the limit of large Re {Nun- <^ Re/Rx] ujy = [{N — l)/N][e'^/CkBT]). The 
result reads: -^ = 1 — ui^g{v) with Gt = 1/NRt, v = eV(3/N and g{x) = [xsinhx — 
4sinh {x/2)\ /8sinh (x/2). Comparison of the above expression with the corresponding one 
derived for the perfectly conducting environment {Re = 0) in [^, -^ = 1 — ^^^u^giv), 
indicates that (^) = t7^(4^) , where -^ = 1 — ^^ stands for the depth of the 

conductance dip. Here, again, the half-width has a universal value given by V1/2 = V^i/2,o- 
The above equation for the asymptotic values of AG/Gt together with the last expres- 
sion for V^i/2 confirms that the effect of a dissipative environment on conductance becomes 
increasingly suppressed by large N and is most noticeable for a two-junction array. 

Finally, let us have a closer look at the results above in view of Coulomb blockade ther- 
mometry (CBT) which is a primary (and secondary) thermometer based on single electron 
charging effects in arrays of tunnel junctions 0. The main parameter is the half- width of 
the CC dip. We conclude that the inaccuracy in temperature measurements arising from 
environmental effects can be made small by increasing the number of junctions, A^. Our nu- 
merical simulations along with the experimental results depicted above show that for A^ = 20 
the temperature determined by CBT is very close to the thermodynamic temperature. In 
practice, in a sample with no intentional Re the effective value of the external impedance is 



of the order of free space impedance Re,eff < -^o — 377 Q, and therefore the agreement is 
better than ±0.5% (see the experimental data point marked by an open circle in Fig. 2). 

In summary, we have studied the effect of the resistive electromagnetic environment on 
transport in arrays of normal metal tunnel junctions within the high temperature and the 
weak tunneling regime. Special attention has been paid to the half-width of the conductance 
curve. Overall agreement between numerical results based on the extension of the PC theory 
and experimental data has been observed. We cannot explain the quantitative discrepancy 
between theory and experiment for the strong enhancement of the width at intermediate 
values of Re in the N = 2 case: higher order tunnelling effects have not, however, been 
included in the present theory. As a practical conclusion, we verify that the effect of envi- 
ronment is mostly emphasized in a two-junction array and can be made sufficiently small 
for thermometric applications by increasing the number of junctions in the array. 

We thank A. Korotkov and J. Konig for discussions. We thank the National Graduate 
School in Materials Physics for support. 



TABLE CAPTION Table I. (a) Parameters of the measured samples with N = 2. 
Samples with the same capital letter, e.g., 2 A, 8A, 9 A and lOA, belong to the same chip, 
(b) and (c) Parameters of the measured samples with A = 8 and A = 20, respectively. 
Samples 4C and IC, 7C and 2C, and 5B and IB constitute pairs of samples, with a space 
of only about 3 fim. between the respective arrays. Capacitances are obtained by fitting the 
theoretical conductance curves to the experimental ones. The only fitting parameter is the 
capacitance of the sample. 



FIGURE CAPTIONS 

Fig. 1. An array of N tunnel junctions in an electromagnetic environment together with 
its bias circuitry. For a symmetric array in a purely dissipative environment, Zej{uj) = Rj, 
and with negligible stray capacitances, Co,k = 0, the real part of the total impedance as seen 
by the A;th junction, Re[Z^^(c(j)], reduces to the simple form presented in the text (Eq. (6)). 

Fig. 2. The measured conductance curve for a two-junction array with Re = 1.12 kQ 
(sample 9A). The curve under the data points is the theoretical conductance curve with 
C = 2.25 fF (see text). 

Fig. 3. Measured half-width of the conductance curve, normalized by bASdNksT/e (see 
text), for samples with N = 2 (solid diamonds), A = 8 (open triangles) and A = 20 (open 
squares) as a function of extra (i.e., that in addition to the unintentional space impedance of 
~ 100 Q) external resistance. For comparison, a sample (A = 20) without any intentional 
on-chip impedance has been shown in the plot by an open circle. The lowermost solid curve 
is obtained by simulation for completely symmetric 20-junction array with C = Ck = 5.0 
fF. The uppermost solid curve is the result of simulation for a two-junction array with 
C = 2.2 fF, whereas dotted and dashed lines correspond to C = 5.0 fF and C = 1.4 fF, 
respectively. The middle solid curve indicates the result of simulation for an eight-junction 
array with C = 5.0 fF. The capacitances were obtained from the depth of the dip of the 
conductance curves. The inset shows the normalized half-width together with the depth of 
the conductance curve {AG/Gt), obtained from the theory, for a 20-junction array. 
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(a) 



Sample 


Re (kfi) 


Rt (kQ) 


cm 


1 


0.10 


39.4 


1.9 


2A 


0.25 


24.7 


3.8 
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0.26 


26.8 


3.5 


4C 


0.35 


15.8 


5.3 


5B 


0.37 


23.0 


3.9 


6F 


0.66 


35.2 


2.4 


7C 


0.83 


13.4 


5.2 


8A 


0.94 


44.4 


2.9 


9A 


1.12 


44.9 


2.3 


lOA 


1.17 


27.4 


3.0 


IIB 


4.76 


52.0 


2.5 



(b) 



Sample 


Re (kQ) 


Rt (kQ) 


cm 


IB 


0.32 


110.1 


3.3 


2B 


3.11 


268.6 


2.3 


3F 


3.48 


170.8 
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(c) 



Sample 


Re (kQ) 


Rt (kQ) 


cm 


IC 


0.70 


145.6 


5.5 


2C 


1.30 


103.1 
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3C 


2.14 


117.5 


7.1 


4D 


2.31 


125.6 


6.1 


5D 


2.49 


122.9 


5.9 


6E 


2.56 


217.2 


3.4 


7E 


2.91 


193.9 


4.0 
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